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ABSTRACT 

Aims. At bright radio powers (/’i. 4 ghz > 10^^ W/Hz) the space density of the most powerful sources peaks at higher redshift than 
that of their weaker counterparts. This paper establishes whether this luminosity-dependent evolution persists for sources an order of 
magnitude fainter than those previously studied, by measuring the steep-spectrum radio luminosity function (RLE) across the range 
lO^'^ < Pi. 4 GHz < 10^* W/Hz, out to high redshift. 

Methods. A grid-based modelling method is used, in which no assumptions are made about the RLE shape and high-redshift be¬ 
haviour. The inputs to the model are the same as in |Rigby et al.| ( |201B : redshift distributions from radio source samples, together 
with source counts and determinations of the local luminosity function. However, to improve coverage of the radio power vs. redshift 
plane at the lowest radio powers, a new faint radio sample is introduced. This covers 0.8 sq. deg., in the Subaru/XMM-Newton Deep 
Field, to a 1.4 GHz flux density limit of S i. 4 ghz ^ 100 ply, with 99% redshift completeness. 

Results. The modelling results show that the previously seen high-redshift declines in space density persist to Ri 4 ghz < 10^^ W/Hz. 
At Ri. 4 ghz > 10“® W/Hz the redshift of the peak space density increases with luminosity, whilst at lower radio luminosities the position 
of the peak remains constant within the uncertainties. This ‘cosmic downsizing’ behaviour is found to be similar to that seen at optical 
wavelengths for quasars, and is interpreted as representing the transition from radiatively efficient to inefficient accretion modes in the 
steep-spectrum population. This conclusion is supported by constructing simple models for the space density evolution of these two 
different radio galaxy classes; these are able to successfully reproduce the observed variation in peak redshift. 
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1. Introduction 

Radio-loud active galactic nuclei are a key component driving 
galaxy evolution; the feedback their expanding radio jets pro¬ 
vide is essential for preventing large-scale cluster cooling flows 
and halting the growth of massive elliptical galaxies (e.g.|Fabian 


et al.|2006l|Best et al.|2006[|Best et al.|2007l[Croton et aL|2006t 


Bower et al. 20061. To understand the timescales upon which 
these processes occur, it is important to first understand the evo¬ 
lution of the radio luminosity function (RLE) to high-redshift. 
An early measurement of this came from [Dunlop & Peacock] 
( |1990| l, who found increases in the space density of both flat 
and steep-spectrum radio AGN of two to three orders of mag¬ 
nitude, over that seen locally. They also saw the first indication 
of an expected higher redshift density decline at z ~ 2.5, corre¬ 
sponding to the build-up of these objects in the early Universe. 
However, their work, and that of subsequent studies (e.g. averj 
et al. 1996[|Jm'vis et al. 200 1[ Waddington et al. 200 l| l, lacked the 
necessary depth and volume needed to unambiguously measure 
this high-redshift behaviour. 

This situation improved with the development of the Com¬ 


bined EIS-NVSS Survey of Radio Sources (CENSORS; Best et 
al.|2003): a survey designed to maximise the coverage of steep- 


RLE. [Rigby eFaL] ( (20TT| hereafter Rll) used CENSORS, com¬ 
bined with additional radio source samples, source counts and 
determinations of the local RLE, to investigate the space den¬ 
sity evolution of the Ei. 4 ghz > 10^^ W/Hz steep-spectrum pop¬ 
ulation via grid-based modelling with no prior assumptions in¬ 
cluded about the high-redshift behaviour. This robustly identi¬ 
fied the post-peak space density decline in the RLE, and found 
that this turnover appears to be luminosity-dependent; at lower 
radio powers (Pi. 4 ghz = W/Hz) the space densities peak 

at z > 1, but the peak moves to higher redshift for the more lu¬ 
minous objects (z > 3 for Ei. 4 ghz > 10^^ W/Hz). A luminosity 
dependence in the position of the steep-spectrum RLE peak can 
be interpreted as a sign of ‘cosmic downsizing’, in which the 
most massive black holes form at earlier epochs than their less 
massive counterparts. This has also been seen for other AGN 
populations, selected at other radio, optical, far-infrared and X- 
ray wavelengths (e.g. De Zotti et al.|2010[[Hasinger et al.|2005[ 


Richards et al.| 2005 [ [Me Alpine et al. 2013| 
2014|l, as well as reproduced in simulations of 


Delvecchio et al. 


rlack hole growth 


(e.g. Fanidakis et al.|20T^ Hirschmann et al.|2012 2014| l. 


spectrum radio sources close to the high-redshift break in the 


Steep-spectrum radio sources can be split into two distinct 
populations: typically powerful objects with ‘standard’ accretion 
of cold gas and likely to be merger driven (‘cold-mode’); and 
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typically weaker sources with radiatively inefficient accretion di¬ 
rectly from their hot gas haloes (‘jet-mode’ or ‘hot-mode’). Lo¬ 
cally, the latter dominate below Pi. 4 ghz ~ 10^® W/Hz, but both 
classes are present at all powers ( |Best & Heckman]|2012[ ). The 
major limitation of R11 was the lack of constraint at lower radio 
powers {P < 10^® W/Hz) out to z > 2. The R11 power range was 
thus only able to probe the RLF regimes where the cold-mode 
sources are in the majority, and therefore could not draw any 
conclusions about the relative contributions of the two classes 
across cosmic time. This paper builds and expands on this ini¬ 
tial foundation by including radio powers an order of magnitude 
fainter than Rll, allowing the RLF evolution to be measured 
across 24 < logPi 4 GHz ^ 28, to high redshift. This improved 
coverage means that the apparent luminosity-dependence of the 
peak space density can be investigated for the first time from 
the brightest sources, down to the regime where the dominant 
steep-spectrum population is shifting between the two types. 

The layout of the paper is as follows: Section describes 
the modelling method used to determine the best-fitting steep 
spectrum RLFs; Section|^presents the SXDF radio source sam¬ 
ple; the results of incorporating this into the modelling pro¬ 
cess are given in Section]^ Section [^discusses the luminosity- 
dependence of the position of the RLF peak; finally the re¬ 
sults are summarised and interpreted in Section Through¬ 
out this paper values for the cosmological parameters of Flo = 
70 km s 'Mpc ', Qm = 0.3 and = 0.7 are used and the radio 
spectral index, a is defined as Sy v““. 

2. Modelling process 

The modelling method in Rll makes no assumptions about the 
shape and evolution of the RLF. It takes as input five radio galaxy 
redshift distributions, measurements of the radio source counts 
over 0.05 mJy < S\aghz 94 Jy, and determinations of the local 
radio luminosity function. The radio samples are selected such 
that they provide a good coverage of the radio luminosity - red¬ 
shift plane out to z ~ 4 (Figure[T]). 

These input data are used to constrain co-moving space den¬ 
sities, p, determined at various points on a P-z grid of radio 
luminosities and redshifts, from which the cosmic evolution of 
the radio galaxy population can simply be read directly. This 
grid is made up of three components - steep-spectrum, flat- 
spectrum and star-forming sources - representing the different 
radio source populations present. Since the steep-spectrum com¬ 
ponent is the dominant population at the luminosities and red- 
shifts considered in Rll, it is the only one allowed to vary in 
the modelling process (see Section[3 for further discussion). The 
constant star-forming grid is created by evolving the local star¬ 
forming galaxy luminosity function of [Sadler et al.|p002| l by 
Pjil + z)^® to z = 2 and by 0.067* at z > 2; the constant 
flat-spectrum grid is taken as the median of the set of evolu¬ 
tionary models presented in [Dunlop & Peacock] ( |1990| l for this 
source population. The starting point for the steep-spectrum grid 
is formed by evolving the local AGN RLF ( [Sadler et al.|2002') l 
by (1 + z)^ in density. 

The grids cover a range of 19.25 < logRi^GHz ^ 29.25 
in radio luminosity, equally spread in steps of AlogPi 4 GHz = 
0.5, and 0.1 < z < 6.0 in redshift, evaluated at z = 
0.1,0.25,0.5,1.0,2.0,3.0,4.0,6.0. The space density at any in¬ 
termediate P-z grid point can simply be interpolated from its 4 
neighbours. Figure [T] illustrates the distribution of grid points, 
together with the corresponding constraints offered by the input 
data. Unconstrained densities are excluded from the fitting pro¬ 
cess. 
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Fig. 1. The distribution of grid points in the P-z plane with the con¬ 
straints offered by the three different types of input dataset highlighted. 
Grey dots show the positions of the individual sources in the four origi¬ 
nal radio source samples taken from Rll. The sources in the new SXDF 
sample are highlighted with red crosses. 



The best-fitting space density grid is determined using the 
amoeba algorithm for downhill simplex minimisation ( [Nelder] 
& Mead|[l965[ ), run in a multi-stage loop with varying scaling 
and tolerance parameters for the initial steps. The goodness of 
fit for a particular model grid is the combined likelihood of the 
predicted source counts, redshift samples and local luminosity 
functions, compared to the real data. The uncertainties in the 
model grid are taken as the marginalized error, calculated using 
the inverse Hessian matrix. 

Full details of both the modelling process and the input 
datasets used can be found in R11. 


3. New input data: the Subaru/XMM-Newton Deep 
Field radio sample 

The faintest radio galaxy sample used in R11 was the AGN sub¬ 
sample of the VLA-COSMOS survey ( [Smolcic et al.|2008] l, with 
a 1.4 GHz flux density limit of 100 //Jy. This dataset was only 
complete to z < 1.3 which, since the next faintest sample had 
a 2 mJy flux density limit, severely limited the high-redshift 
constraint in the P-z grid at low radio powers. This, therefore, 
restricted the RLF evolution analysis to logPi 4 GHz ^ 25 only. 
To improve this constraint a radio survey of comparable depth 
but better redshift coverage, carried out in the Subaru/XMM- 
Newton Deep Field (SXDF; Simpson et al.|20()6 |2012| l, is now 
included in the modelling. This covers 0.81 sq. deg. and contains 
a complete sample of 505 sources with S i. 4 ghz ^ 100 //Jy. Ro¬ 
bust redshifts, both spectroscopic (51%) and photometric, exist 
for 99% of the sample. The photometric redshifts for the remain¬ 
ing 7 sources are contaminated by foreground objects and are 
therefore excluded here. 

At this flux density level >50% of the sources are expected to 
be star-forming galaxies (e.g. Norris et al.|2013 1 . Since this pa¬ 
per is concerned with the behaviour of the steep-spectrum pop¬ 
ulation (which comprises the majority of the remainder), the cat¬ 
alogue is restricted to radio-loud AGN only, as described below. 
In principle, this step is not necessary as the evolution of the star¬ 
forming RLF in the modelling process should be accounted for 
by the treatment applied to the star-forming grid. However, un¬ 
certainties in the evolution parameters used for this, combined 
with the non-negligible contribution of this population in this 
survey, means that it is prudent to perform this step to limit the 
effect of any errors. 
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Fig. 2. The redshift distribution histogram for the SXDF catalogue after 
removing star-forming galaxies. The solid line shows the best-fit model 
prediction for this sample, and the dot-dash line shows the 12th order 
polynomial used to represent the distribution in the modelling process. 
The contribution to the best-fit by the flat-spectrum grid is also shown 
(dotted line) to illustrate the dominance of the steep-spectrum grid at 
all redshifts. 


The AGN/star-formation separation is done by calculating 
the ratio of mid-infrared, S 24 tim, to radio flux density via the k- 
corrected q parameter: q - \Qg{S 2 AiimlS i. 4 ghz); following |Don- 
|ley et al.] ( |2005| l, AGN are defined as sources with q < 0, in- 
dicating excess radio emission over that from star-formation. 
Applying this cut results in a final sample of 179 sources. The 
q parameter adopted is higher than the q - -0.23 value used 


previously to select radio-loud objects (e.g. Ibar et ah] 
However, inspection of Figure 6 in |Simpson et al.| ( |2012^ shows 
that the higher cut is more appropriate for this sample. In prac¬ 
tice, the difference in the modelling results from changing q is 
negligible within the errors. 

The AGN cut also removes radio quiet quasars from the sam¬ 
ple. These also begin to be a larger fraction of the radio popula¬ 
tion at S 1 . 4 GHz < 1 mJy, and display radio emission which may 
arise mainly from star-formation instead of black hole accretion, 
though this is still a topic of active debate ( [Smolcic et al.|2015| 
and references therein). Their contribution to the overall popu¬ 
lation is not explicitly included here, but rather assumed to be 
covered by the star-forming grid. 

The redshift distribution of the sample is shown in Figure]^ 
together with the 12th order polynomial used to represent it in 
the modelling process. This step is necessary as the low number 
of sources at high-redshift makes a direct calculation of the 
difficult (see R11 for further details). 

The VLA-COSMOS sample previously used in the mod¬ 
elling process is now replaced here by the SXDF sample de¬ 
scribed above. 


4. Model results 

4.1. Dataset comparison 

The independent redshift distribution generated from the best-fit 
steep-spectrum P-z grid is shown in Figure this agrees well 
with the real dataset at z > 0.3. Below this the SXDF potentially 
misses resolved sources ( Simpson et al.|[2012[ ), which may ex¬ 
plain the apparent overprediction of the model at low redshift. 


The Figure also shows the small relative contribution to the total 
population from the un-fitted flat-spectrum grid. 

The total number of AGN (i.e. flat- plus steep-spectrum 
sources) given by the model is 199 ± 31 (188 + 29 at z > 0.3), 
which compares well to the actual figure of 179 (178 at z > 0.3). 
In addition, the predicted number of star-forming galaxies (377) 
broadly agrees with the 326 observed in the SXDF sample and 
suggests that the treatment of this population in the modelling 
is sufficient to not give rise to additional uncertainty in the grid 
results. 

The agreement between the model and the other input 
datasets is also good, however these comparisons are not ex¬ 
plicitly shown here as they are essentially unchanged from the 
results previously presented in R11. 

4.2. Model RLFs 

The predicted variation of space density with redshift for 8 radio 
power bins (24 < log P < 28) is shown in Figure|^ Note that the 
p for a particular P-z point is only shown if it is constrained by 
both the source counts and the redshift distributions (Figure[2l. 

The main benefit of including the SXDF sample in the mod¬ 
elling comes at the lowest powers (logP < 25.5), where it al¬ 
lows the RLF behaviour to be studied at z > 1. Inspection of 
Figure shows that decreases in space density are also present 
here, though these tend to occur at lower redshifts than for higher 
power ranges; this will be discussed further in Section]^ 

As part of the modelling process the original P-z grid is con¬ 
verted into a grid of flux density and redshift to allow easier 
comparison with the source counts and redshift distributions. To 
make this 5-z grid a spectral index, a, is needed; following R11 
this is assumed to vary with redshift as cr = 0.83 + 0.4 log(l -H z) 
( |Ubachukwu et aL]|1996| l, with an additional uncertainty of 
a + 0.2 at each redshift to account for any inter-bin variations. 
Figure shows the range of best-fit results returned when this 
assumption is changed to either a constant, but very steep spec¬ 
trum {a = 1.5), a constant value of 0.8 (the mean value at low 
redshift), or a stronger increase with redshift {a = 0.8 -H 0.25z). 
The general effect of using the steeper a values is to increase 
the overall density values, and weaken, though not remove, the 
z > 1 space density turnovers as sources move to bins in the S -z 
grid corresponding to higher radio powers. The assumption of 
the shallower value of a = 0.8 has the opposite effect. The small 
spread in results seen at 26 < logP < 26.5 following this a 
variation suggests that this is the radio power bin which is most 
strongly constrained in the modelling process. 


5. The luminosity-dependence of the high-redshift 
turnover 

Inspection of Figure|^suggests that, as in R11, the position of the 
space density peak is luminosity dependent. This can be clearly 
seen if the model steep-spectrum RLFs are shown together, as 
in Figure]^ The addition of the SXDF sample as a constraint to 
the RLF modelling allows this behaviour to be probed to lower 
radio powers than previously possible. 

The position of the peak redshift in each radio power range 
is determined by fitting a polynomial (generally of order 2, but 
of order 3 where necessary) to the model space density distri¬ 
butions in the P-z grid. This overcomes the ambiguity arising 
from the discrete nature of the P-z grid, and the wide redshift 
bins used. Figure]^ shows the results of this for the best-fitting 
steep-spectrum grid presented in Section |4.2| together with the 
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24.00 < log P < 24.50 




0.1 1.0 10.0 


27.00 < log P < 27.50 



24.50 < log P < 25.00 



25.50 < log P < 26.00 



26.50 < log P < 27.00 



27.50 < log P < 28.00 



Fig. 3. The model steep-spectrum RTFs, as a function of redshift, from the best-fitting P-z grid (black shaded region). Results are only plotted if 
they are constrained by at least two of the input datasets (Figure[^. Dark grey shaded regions show the spread in results from varying the assumed 
value of the spectral index. Dotted and dashed lines show the predictions of the simple RTF model described in Section [F!2| for the hot-mode and 
cold-mode populations respectively (for clarity these are plotted up to the peak redshift value only). The brightest RTFs (log P > 26) were fully 
constrained previously, and are therefore unchanged from R11. They are reproduced here for completeness and comparison with the simple RTF 
model. 


spread arising from repeating this process for the grids deter¬ 
mined using altered spectral indices. The small variation in val¬ 
ues at logP = 26.25 supports the previous identification of 
26 < logP < 26.5 as the most constrained radio power bin in 
the model grid. 


The luminosity-dependence of the peak redshift, Zpeak, is 
clear, with the most powerful radio sources peaking at earlier 
times than their weaker counterparts. However, the extension to 
lower radio powers reveals an apparent flattening of the rela¬ 
tionship; below logP < 26.5 the results are consistent with a 
constant value of Zpeak- 
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Fig. 4. The model steep-spectrum RLFs, as a function of redshift, from 
Figurej^plotted together to illustrate the changing position of the space 
density peak. Larger bins of A log P = 1 are used for clarity. 
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Fig. 5. The variation in the redshift of the peak space density with radio 
power for the best-fitting steep spectrum grid determined here (black 
line), compared with that predicted using the simple model described in 
Section |5^ (blue line). Dotted lines delineate the relative contributions 
to the latter from the two underlying populations. The error bars show 
the uncertainty in the polynomial fits used to determine the peak posi¬ 
tion, whilst the shaded region represents the range in results found from 
varying the input parameters used in the RLF grid modelling. 


5.1. Comparison with quasar behaviour 


The luminosity-dependent behaviour of the RLF peak redshift 
seen here for steep-spectrum radio sources has also been ob¬ 
served in other AGN populations (e.g. De Zotti et al. 2010t 
Hasinger et al.|[2()05[ |Richards et al.|[2()05t |Delvecchic> et al.| 
2014|l. A direct comparison of the Zpeak changes is possible 
for tlat-spectrum quasars, using the luminosity function (QLF) 
models of [Hopkins et al.| ( |2007| l; these are derived from a large 
set of observed QLFs, from the mid-infrared to X-ray. 

The optical B-band (0.44 pm) is chosen for extracting the 
quasar behaviour, for ease of comparison to the luminosities in 
the radio waveband. It should be noted however that the variation 
in peak redshift calculated using the bolometric QLF instead is 
qualitatively similar. The Hopkins et al. model is queried for a 
luminosity range over which it is reliable (20 < log Pu-band < 
27), and the Zpeak value determined in each case. The subsequent 
result is shown in Figure together with that already found for 
the steep-spectrum sources. The ratio of the two curves is also 
shown, as a function of peak redshift. 



^peak 


Fig. 6. A comparison of the variation of the redshift of the space den¬ 
sity peak with radio luminosity (steep-spectrum sources; black solid 
line, with the uncertainty range arising from the assumed spectral in¬ 
dex indicated by dotted lines) and B-band luminosity (quasars; blue 
dashed line). The subplot shows the ratio between the two luminosities 
as a function of Zpeak- The dotted line here indicates a constant ratio of 
log(Pi . 4 GHz/PB-band) = 1-5 - the minimum value historically adopted 
when classifying a source as radio-loud ([Kellermann et al.|1989^. 


The shape of the Zpeak-log PB-bani/ curve derived for the 
quasar population is broadly similar to that of the radio galaxies, 
with indications of flattening at 24 < log Py < 26, before con¬ 
tinuing to decrease at lower powers. The log(Pi. 4 GHz/PB-band) 
ratio begins to increase from a constant value of ~ 1.5 at 
z > 3, within the uncertainties; it is interesting to note that 
log(Pi . 4 GHz/PB-band) ^ 1.5 is the historical definition of a radio- 
loud quasar (Kellermann et al.|1989|l. 


5.2. The reiation to cosmic downsizing 


The obvious physical interpretation of the luminosity-dependent 
space density peak seen here is cosmic downsizing, in which the 
most massive black holes have formed by z ~ 4, with the less 
massive forming at later times. This apparently contradictory 
anti-hierarchical growth can be explained if the dominant mech¬ 
anism of AGN fuelling changes with cosmic time from cold gas 
accretion via major mergers (‘cold-mode’) to radiatively ineffi¬ 
cient accretion directly from hot gas haloes (‘jet-mode’ or ‘hot¬ 
mode’). Simulations of AGN evolution support this picture; they 
show that hierarchical black hole growth is able to produce the 
observed downsizing trend if inefficient accretors come to dom¬ 
inate at low redshift (jFanidakis et al.||2012[ [Hirschmann et al.j 
[20T4l l. 


Locally, hot-mode sources make up the majority of the 
radio-galaxy population at logPi 4 GHz ^ 26 (Best & Heckman 
|2012| l, though both types are generally found over the range of 
radio powers studied. In comparison, by z ~ 0.7 the space densi¬ 
ties of the two are similar at low radio powers, with the cold¬ 
mode sources becoming more dominant at logP 


1.4GHz 


> 26 

[Best et al.||2014|l. The flattening seen here in the Zpeak vs. 
ogPi 4 GHz relation occurs at corresponding radio power to this 
transition, lending further support to the assertion that the RLF 
evolution seen here directly represents the interplay between 
the two populations. It also suggests the weaker, likely hot¬ 
mode, sources undergo little evolution in space density; a finding 
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Table 1. The peak redshift (Zpcak) values used in the simple model of hot and cold-mode RLF evolution, described in Section 


5.2 


In this model 


hot-mode sources increase as (1 +z) to Zpeak = 0.8, whereas cold-mode sources have a luminosity dependent value of Zp^ak and a pre-peak increase 
of (1 + z)^. Both populations have a post-peak decline of (1 + The numbers in brackets are the peak values after the predicted RTFs are 

re-sampled using the same redshift bins as in the grid modelling. The ‘Hot+Cold Mode’ column gives Zp^ak for the combined hot and cold-mode 
RLFs, extracted via polynomial fitting. The ‘Grid Model’ column gives the corresponding Zpcak determined from the best-fit model grid (Section 
|4.2j >; the lower and upper bounds give the range arising from varying the input parameters used in the grid modelling, as shown by the shaded 
region in Figure]^ 


Radio power range 



^peak 



Hot-mode 

Cold-mode 

Hot-bCold mode 

Grid model 

24.00 - 24.50 

0.8 (0.5) 

1.8 ( 1 . 0 ) 

0.4 

1 # 

Q-3 

24.50-25.00 

0.8 (0.5) 

2.0 ( 2 . 0 ) 

1.1 

25.00-25.50 

0.8 (0.5) 

2.2 ( 2 . 0 ) 

1.2 

25.50-26.00 

0.8 (0.5) 

2.3 (2.0) 

1.3 

1 2 ”-^ 

Q-3 

26.00 - 26.50 

0.8 (0.5) 

2.4 (2.0) 

1.3 

1 40.0 

Q -8 

2 3°-^ 

Q-4 

3 4° ' 

_-’•^ 1.2 _ 

26.50 - 27.00 

0.8 (0.5) 

2.7 (2.0) 

1.3 

27.00 - 27.50 

0.8 (0.5) 

3.0 (3.0) 

2.2 

27.50-28.00 

0.8 (0.5) 

3.6 (3.0) 

2.2 


broadly replicated in Best et all] ( 2014| l which directly catalogued 
sources into the two classes. 

This assertion can be tested by predicting the dual popula¬ 
tion space density evolution, taking as a starting point the well- 
defined local RLFs for the two classes ( [Best & Heckman|2012|l. 
For the evolution of the hot-mode population, we follow |Best^ 
al. (2014|l in assuming that their evolution broadly follows that 


of the space density of massive quiescent galaxies, which consti¬ 
tute their host galaxies. [Best et aL] ( |2014) l showed that the space 
density of these galaxies remains relatively flat out to z ~ 0 . 8 , 
and then rapidly decreases as (1 We modify this to allow 

for an increase in the space density as (1 H- z) out to the peak 
redshift at z = 0 . 8 , as this is a better match to the evolution of 
the 10 ^“* < Pi. 4 GHz < 10^^ RLFs in the model grid. 

For the cold-mode population, we determine the 
(luminosity-dependent) peak redshift by comparison with 
the optical quasars. Figure suggests that there is a roughly 
constant ratio between the luminosities of radio galaxies and 
quasars at high values of Zpeak, where the cold-mode sources 
dominate. Assuming this is true for this population, Zpeak for a 
particular radio luminosity is therefore taken as the peak of the 
corresponding P-band quasar luminosity. For this population, 
we allow a more rapid pre-peak increase in space density of 
(1 H- z)^, which matches the evolution seen at the highest radio 
luminosities in the model grid. The post-peak decrease is again 
taken as (1 -H z)“®'^. 

The predicted RLFs for the two populations are then com¬ 
bined, sampled using the same luminosity and redshift bins as in 
the model grid, and the peak redshift extracted. 

The variation in peak redshift with radio luminosity arising 
from this simple model is shown in Figure It does well at re¬ 
producing the overall behaviour derived from the best-fit model 
grid. At the faintest (logPi^oHz < 24.5) luminosities the pre¬ 
dicted RLFs are dominated by the hot-mode sources, whilst the 
cold-mode sources dominate at logPj 4 ghz > 26.5. In-between, 
the RLFs are a combination of the two classes. This is reflected 
in the changing shape of the Zpeak - log F’i. 4 ghz curve; the plateau 
at 24.5 < logf’i. 4 GHz < 26.5 arises, in this scenario, from the 
transition from the cold-mode to the hot-mode populations. 


6. Summary & Conclusions 

The results presented in this paper have demonstrated the ben¬ 
efits of extending the P-z plane coverage to lower radio pow¬ 
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ers when measuring the evolution of the steep-spectrum RLF. 
The addition of the updated SXDF sample to the constraining 
datasets used in the modelling process shows that the high- 
redshift turnover in the RLF persists to logPi 4 GHz > 24 - an 
order of magnitude fainter than could be measured previously. 

The redshift at which the peak space density occurs is found 
to be luminosity-dependent at logPi 4 GHz ^ 26, with the most 
powerful sources peaking at earlier times than their weaker 
counterparts. This mirrors that seen when considering quasar op¬ 
tical luminosity functions, which is to be expected as the radio 
population is dominated by cold-mode (i.e. quasar-like) sources 
at these luminosities. The presence of a mean ratio for these pow¬ 
ers (log[Pi 4 GHz/^B-band] ~ 1-5) which remains constant with 
Zpeak again suggests that the radio and optical observations are 
seeing the same objects. This picture changes at lower radio lu¬ 
minosities (logPi. 4 GHz ^ 26). Flere the value of the peak redshift 
remains constant within the uncertainties, whereas it continues 
to fall for the optical population. This implies an increase in the 
F’i. 4 GHz/F’B-band ratio and suggests that the dominant radio pop¬ 
ulation is changing from cold-mode to hot-mode. 

Simple models for the space density evolution of the hot¬ 
mode and cold-mode populations are able to reproduce the ob¬ 
served luminosity dependence in the peak redshift, reinforcing 
the conclusion that this behaviour arises from the transition be¬ 
tween the two populations. To properly understand the relative 
contributions of the hot-mode and cold-mode sources to the 
RLF evolution, their behaviour should be minimised separately 
in the fitting. This is trivial to do with this grid-based modelling 
method, but the input data needed are currently lacking across 
the required redshift range. Future large, deep, radio surveys in 
regions with extensive multiwavelength coverage will make this 
possible, whilst also allowing the coverage of the P-z plane to 
extend to radio powers of logPi. 4 GHz < 24. 
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